RKKY interaction in the nearly-nested Fermi liquid 
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We present the results of analytical evaluation of the indirect RKKY interaction in a layered metal 
with nearly nested (almost squared) Fermi surface. The final expressions are obtained in closed 
form as a combination of Bessel functions. We discuss the notion of the "2A;f" oscillations and show 
that they occur as the far asymptote of our expressions. We show the existence of the intermediate 
asymptote of the interaction which is of the sign-reversal antiferromagnetic type and is the only term 
surviving in the limit of exact nesting. A good accordance of our analytical formulas with numerical 
findings is demonstrated until the interatomic distances. The obtained expressions for the Green's 
functions extend the previous analytical results into the region of intermediate distances as well. 



75.30. Et, 74.72.-h, 75.20.Hr 

The Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action was found to play an important role in various 
problems involving the interaction of localized moments 
in a metal via polarization of conduction electrons. The 
spatial dependence of this interaction for the spherical 
Fermi surface (FS) in three dimensions was obtained in 
Rcf.tl about 40 years ago. It was demostrated there that 
at large distances r the interaction decays as r~ 3 and has 
the 2kp oscillations withXhe Fermi momentum kp. Later 
Roth, Zeiger and Kaplana have generalized this result for 
the case of non-spherical FS. The interaction was repre- 
sented as a series in 1 Jr and the existence of the direction- 
dependent period of oscillations (2k' F )~ 1 was shown. 

The limitation of this latter result is the inapplica- 
bility of the theory at small distances k* F r < 1, i.e. in 
the region where the RKKY interaction is mostly signif- 
icant. It could be among the plausible reasons for the 
nowaday situation, when the notion of 2kp oscillations 
is widely explored while the actual meaning of this term 
for the non-spherical FS remains uncleat,in many cases. 
Numerical calculations are instrumentaH to extend our 
understanding of non-spherical FS ever since 1957, but 
a theoretical understanding even on a qualitative level is 
very important. 

In a recent paper one of the authors presented the 
method which enabled to obtain in a simple manner 
the RKKY interaction for the spherical FS and any di- 
mensionality of the system. The closed expressions were 
found both in r— and q— representations^ 

In the present paper we generalize this method to ar- 
rive to the closed analytical expressions for the RKKY 
interaction in a layered (2D) metal with highly non- 
spherical (nearly nested) Fermi surface. Our main results 
could be summarized as follows. 

i) The RKKY interaction for this type of the FS consists 
of three parts. The first one stems from the flat regions 
in the fermionic spectrum. This term decays mostly as 
1 Jr and is important only along the diagonals x — ±y in 
the r— space. 

ii) The second part of the RKKY interaction is con- 
tributed by the vicinities of the points (0, tt) and (n, 0). 



These are the saddle points in the fermionic dispersion 
giving rise to the van Hove singularities. At r — > oo this 
part of RKKY takes its exact correspondence with the 
previous findings .□ It turns out however that it is the far 
asymptote of the interaction. 

iii) At k F r < 1 the third part of the interaction comes 
into play. This intermediate asymptote arises as an in- 
terference between contributions from the vicinities of 
the different van Hove points. This term possesses the 
overall prefactor cos(Qr) with the spanning wave-vector 
Q = (tt, it). In the nearly nested situation this last term 
is present down to the interatomic distances and strongly 
favors the commensurate antiferromagnetic ordering of 
the localized moments. 

iv) The last but not least. We show both analytically 
and numerically that our expressions are valid near the 
interatomic distances. This result could be qualitatively 
explained in the following way. It is clear that to eval- 
uate the interaction at the distances r one should know 
the details of the fermionic dispersion on a scale 1/r in 
k— space. Hence the finest details of the Fermi surface 
are of importance at the largest distances. One may also 
conclude that even rough knowledge of the spectrum over 
the whole Brillouin zone is enough for a good descrip- 
tion of the corresponding quantities at the interatomic 
distances in r— space. This is exactly what our method 
does by grasping the key features of the FS. 

The rest of the paper is organized as follows. We for- 
mulate the problem and introduce the basic ingredients 
of our treatment in Sections I and II. The Green's func- 
tions resulted from the van Hove points of the spectrum 
are found in Section III. Their contribution to the RKKY 
interaction is analyzed in Section IV. The role of the flat 
parts of the spectrum is discussed in Section V. We make 
the concluding remarks in the Section VI of the paper. 



I. GENERAL FORMALISM 

We begin with conventional form of the exchange in- 
teraction between the localized moment J and electron 
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spin density s(r) : 



V(r) = -AJ(R)s(r)<S(R - r) 



(1) 



Here A is the exchange coupling constant. The RKKY 
interaction between two localized moments via the con- 
duction electrons may then be written in the following 
form 



H 



1 



RKKY 



A 2 3 1 J 2X (r ll r 2 ) 



(2) 



where r— dependent part of the interaction coincides with 
the Fourier transform of the non-uniform static suscepti- 
bility x(q) (Lindhard function) and is usually written in 
the form : 



X(ri,r 2 ) 



(2-kY 



d 3 kd 3 q 



^ e i(q-k)(n-r 2 ) 



(3) 



with the unit cell volume vq and the Fermi function 
iik = (exp(efe/T) + 1) . For our purpose, it is more 
convenient however to represent the above expression in 
the equivalent form 



X(ri,r 2 ) 



-T^G(ito n j 1*2, T\ 



(4) 



here Matsubara frequency ui n — irT(2n + 1) and G is 
the electronic Green's function. In the case of low tem- 
peratures considered below we use the limiting relation 

te„-/_°°>/(h 

The electronic Green's function is given by 



(27T) 



d 3 k 



exp(zk(ri - r 2 )) 



U k (Ti)u%(T 2 ) 

(5) 



with the periodic Bloch function Uk (r) . Generally speak- 
ing the u— functions should be inserted into (^), and it 
is the case of free electrons only when the RKKY in- 
teraction depends on the absolute value of the distance 
ri — r 2 . We consider below the tight-binding form of the 
electronic Hamiltonian which implies the following form 
of u k (r]B: 



Ufe ( r ) = H e 



V(r- a„), 



(6) 



here the Wannier function ip(r) rapidly decays with dis- 
tance and is close to the atomic wave- function at small r. 
In view of this rapid decrease on the scale of interatomic 
distances, it is possible to neglect the dependence of Ufc(r) 
on k within the first Brillouin zone for most positions of 
r in the unit cell. Therefore we may replace Ufe(r) by 
uq(t) and assume the exchange coupling is appropriately 
redefined, A — » A|-u ( r )| 2 - It is clear that upon this re- 
definition the Green's function (||) depends only on the 
difference r = ri — r 2 . 



In view of forthcoming consideration we stress that 
the above replacement of the amplitude of Bloch func- 
tion is the only uncontrolled apaoximation of our treat- 
ment. The relevant calculations! show that ip(r) for the 
nearest neighbors is of order of magnitude smaller than 
99(0). Therefore we expect that omitting the dependence 
of Ufc(r) on k might cause only minor corrections to our 
results. 



II. NEARLY NESTED FERMI SURFACE 

We consider a two-dimensional case of almost nested 
Fermi surface, with the quasiparticle dispersion given by 



£k = — 2i(cos k x + cos k y ) — fi, \fi\ -C t 



(7) 



with t > 0. This tight-binding form of the spectrum 
particularly appearecLin different models related to the 
high-T c phenomenon!] Henceforth we let the lattice pa- 
rameter to be unity. 

We will regard the different parts of the spectrum (0) 
on the different manner. We divide the whole FS onto 
four parts, which are as follows. 

The vicinities of the saddle points of the spectrum. The 
dispersion near the points (0, ±7r) and (±7r, 0) are given 
by 

e(o,=br)+k - t{kl - hi) - f i + 0(tk 4 ), (8) 



£(± w ,o)+k - -t{k 2 x -k 2 y )-^ + 0(tk 4 ), (9) 
respectively. 

The flat parts of the Fermi surface. In the vicinity of 
the pair of wave- vectors ±(7r/2, 7t/2) the dispersion takes 
the fojm similar to the one-dimensional case near half- 
fillineQ: 



£k±(7r/2,7r/2) - ±2t(k x + k y ) - fj, + Q(tk 3 ). 



(10) 



The pair ±(tt/2, 
dispersion law 



-7r/2) is characterized by the similar 



Ek±(7r/2 -tt/2) - ±2t(fcx - ky) -fi+ 0{tk 3 ). (11) 



The expansions (@)- (|Il]) are valid until fc < 1. In its 
turn, it means that dropping the higher terms of the ex- 
pansions, one hopes to obtain the correct form of RKKY 
interaction at r > 1. Below we revise this statement and 
compare our analytic results with the numerical findings. 

We see that the whole vicinity of the FS can naturally 
be divided into parts of the mainly two-dimensional hy- 
perbolical and linear character of dispersion. Therefore 
we can represent the Green's function as a sum of eight 
different contributions which is symbolically written in 
the form 



G(iuj,r) = ^G ko (iw,r) 



(12) 
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where k = ±(0, tt), ±(tt, 0), (±7r/2, ±7r/2) and the par- 
tial Green's function G^iito, r) originates from the inte- 
gration over the part of the Brillouin zone with particular 
form of dispersion law : 



G ko (iw,r) 



(27T) 



d2k exp(z(k + k )r)^ 

fc<l zw _ £ k+k 



It should be noted that this decomposition is a general 
one and may be done in all cases when the FS possesses 
the saddle points and flat parts. Hence our consideration 
of the nearly nested FS may be regarded as a particular 
example of the analysis of the RKKY interaction for the 
strongly non-spherical FS. 



III. THE GREEN'S FUNCTION FROM THE 
SADDLE POINTS IN THE SPECTRUM 

In this Section we consider the contribution to the 
Green's function from the vicinities of the saddle points of 
the spectrum only. The flat parts are analyzed in Section 
[v|. As we will show below the saddle points determine 
the RKKY interaction almost at all directions in r— space 
except the regions near the diagonals x = ±y where the 
contribution from the flat parts of the FS becomes im- 
portant. 

At the first step we use the following auxiliary repre- 
sentation of the Green's function : 



G(iw,r) 



(27T) 



dr /d 2 kexp[ikr + re la [iuj - Ek] 



(14) 



where we put a = sign{uS)'K /2. It is evident that e la = 
i sign(uj) but this formal trick facilitates the following 
analysis of the complex- valued expressions. 

The evaluation of the saddle point contribution to the 
Green's funciton can be made in more general form ap- 
plicable for the parts of the. electronic spectrum near the 
so-called stationary pointsel. By definition these are the 
points ko where £k is well approximated by the expres- 
sion 



Ek +k ~ ^kni 1 k 



(15) 



with the generally anisotropic tensor of masses m and 
some effective chemical potential p. We assume the ap- 
plicability of (|l^) for the wave-vectors k not exceeding 
some scale n, comparable to inverse lattice parameter. 
In other words the significant parts of the FS may be 
approximately mapped by ()l5|). 

In view of (15) the integration over k in ( [l4| ) becomes 
Gaussian one and is easily performed. The only compli-, 
cation here comparing to the case of free electron gascl 
is the finite value of k which restricts the validity of the 



subsequent equations by r > k~ 1 ; we discuss it in more 
detail below. 

Thus we obtain the following expression for the partial 
contribution to Green's function from the vicinity of kg 



G k o(^r) = - e -^Me-or 

ZTT 

dr r ia P - 

— expiree e 

t 2r 



(16) 



where z — p + iu> and p = rffir is the square of distance 
in the metrics defined by the mass tensor. The last in- 
tegral is [expressed via the modified Bessel (Macdonald) 
functionp namely 



Gk (iLo,r) = -e 



, a v\d^m e ^ Ko{VZ ^- p) (17) 



In the equation (JT^) the branch of \/~2zp is chosen 
from the condition of its positive real part. In partic- 
ular case of 2pp > 0, this latter condition means that 
the argument of Macdonald function Kq(^/— 2zp) has a 
discontinuity at uj = 0j3 



K {^2z~p) 




u/fx -> +0 



(18) 



( 1 2) 

where ' '(x) are Hankel functions. Note that in the 
case of spherical Fermi surface the condition 2pp = 



> is fulfilled automatically. For the electronic 



type of dispersion m > and p, > while for the hole- 
like dispersion law one has m < and p < 0. It is clear 
from Eq. ( |l8| ) that for non-spherical Fermi surface the ef- 
fective Fermi momentum is given by \f2pp = k* F r and 
may strongly depend on the direction in real space (cf. 
Eq.© below). 

Let us discuss now the region of applicability of the 
expression (|l7f). The Gaussian integration in (Q) for ej. 
given by (|15|)is justified upon two conditions. First, the 
center of quadratic form in k which is r _1 rm should lie 
within the circle of radius k, where the expansion jl5| ) 
is applicable. Second, one should demand the criterion 
r|Km _1 k\ 3> 1 to be satisfied , to ensure the Gaussian 
value (p^). In the principal axes j of m these criteria can 
be combined as follows : 



t S> max 



(19) 



At this point our analysis is somewhat different for the 
cases k* F r = \/2pp ^> 1 and k* F r < 1. 

In the first case of largest r the principal contribution 
to the integral ( |l6| ) is deliveredby r = tq(1 + 0(1/ y / 2zp)) 
with To = \J p/2z = k F r / (2p)[3. Simple arguments show 
then that for the ellipsoidal FS the criterion (|l^) is always 
fulfilled, coinciding with the obvious demand k* F -C k. 
For the spectrum with the saddle point the condition ( |l9| ) 
may be violated at large distances and near the nodes 



3 



p ~ 0. Let the angle if be measured from the i-axis 
in r-plane. Then k F is given by ( ^6|) and we lose the 
applicability of (JT^) at 



« v/|(^±7r/4| < 



(20) 



i.e. within a narrow sector along the diagonals if r/a 3> 

IVmI » i. 

In the case of fc^r < 1 it is the region |p| < r < |z| 
which is essential in (|l^). Substituting the lower bound- 
ary t = |p| into ( |l9| ) we come to the evident condition 
nr > 1 for the ellipsoidal FS. For the FS with the saddle 
point the criterion (19) is again violated near the node of 
p. For our particular form of the spectrum the Eq. ( |l7| ) is 
applicable outside the band along the diagonals x = ±y 



r\(p ± 7r/4| > k 



(21) 



We depict the above regions of applicability on the Figure 

Let us now consider the particular form of the Green's 
functions arising in the case of the tight-binding spectrum 

©■ 

The integration over k being restricted to the first Bril- 
louin zone leads to the additional phase factor in (|D]) as 
discussed in Appendix |a[ The desired expressions ac- 
quire the form : 

G (0>7r) (iu,,r) = Jji^K (^2Tp)^ sisn ^ (22) 



isign{ui) 
2nt 



a i7r\x\sign(uj) 



(23) 



where p = (x 2 — y 2 )/2t and ^/\det in| is replaced by its 
actual value l/2t. 



IV. EVALUATION OF THE RKKY 
INTERACTION 



A. integer values of r 



We note that if the values of coordinates x,y in (22), 
( p3| ) coincide with the integer numbers of lattice periods 
then one has e^M^snM = e i™ ^ e in\y\sign(u) = e iwy and 
e 2nrx _ e 2nry _ ^ ^ e cons i(jer this simpler case first, 
while the case of non- integer x, y is discussed in the next 
subsection. 

Away from the diagonals x = ±y in r-space, we have 
from (0) the following expression : 



X(r) 



1 



"+* 00 dz 



47T 2 t 2 ./ M - ioo 27ri 



K 2 (^2z~p) + K 2 (^2z~p) 
+2e i ^ + ^ (v /= 2^)#o(v / 2^)l (24) 



Here we redefined the variable of integration lu — > z = 
p + iu>. 

As we discussed above if pp > then the function 
Ko(yJ—2zp) has the discontinuity ([l^) at 2 = p. On the 
contrary the second term in ( p4| ) is continuous function 
of z in this case and the corresponding integral is zero, 
because one can shift the integration contour to z — > +oo 
wherein Kq(z) oc e~ z . 

At pp < the situation is reversed, hence one can 
combine these two cases and castlhe contribution of first 
two terms in d24|) into the formtil 



■$i(k* F r) 



$i(o) = J (o)Fo(o) + Ji(o)ri(o) 



(25a) 
(25b) 



with Besscl functions J n (x) and Y n (x) and the direction- 
dependent value of the effective Fermi momentum 



k F — 



VWpp\ 



p 

— cos(2<^) 



(26) 



where we restored in the rhs the lattice parameter a. 
Note that the expression (|25|) gives the RKKY interaction 
for the cylindrical Fermi surface as well, in which case k* F 
coincides with the conventional Fermi momentum!!. The 
only difference is in the general minus sign resulted from 
the sign-indefinite property of the mass tensor, det(m) — 
-l/(4i 2 )<0. 

The asymptotes of this part of RKKY interaction un- 
der criteria ( pC| ) and (|l|) are as follows 



Xi( r ) 



\fi\ sm(2k F r) 
8ir 2 t 2 (k* F r) 2 ' 

4ir 2 t 2 



k* F r > 1 



lnfc^r, k* F r < 1 



(27a) 
(27b) 



Therefore in the limit of large distances the power-law 
decrease of Xi( r ) is accompanied by oscillations with 
the direction-dependent period 2k F , in accordance with 
usual expectations, while for the small k* F r these oscilla- 
tions are replaced by logarithmic singularity. 

The third term in ( |24|) is the integral of the prod- 
uct of the Green's functions resulting from the regions 
of the FS with the different character of dispersion. As 
a result, one has the additional prefactor of the form 
exp(iQor) = (— \) x+ y (for integer x,y). The appear- 
ing "antiferromagnetic" wave- vector Qo = {it, 7r) merely 
connects two regions in the Brillouin zone, where the dis- 
persion is close to the Fermi level. In this third term of 
Eq. (|24|) the discontinuity at z = p exists for berth signs 
of pp. The integration over z can also be donetil and we 
obtain after some calculations : 



X2{v) = e^^<S> 2 (k F r) 



$ 2 (a) 



47T 2 t 2 

J (a)#i(a)- Ji(a)Jf (a) 



(28a) 
(28b) 



4 



The asymptotes of this expression are 
X 2 (r) = e^M^^^e-^, **r » 1 (29a) 



4tt^ 2 {k* F r) 2 



JQ T_W_ 



AttH 2 (k* F r) 2 ' 



ktr < 1 



(29b) 



Let us discuss the significance of this second part of 
RKKY interaction. First, one sees from (p7|), ( |29| ) that 
in the far asymptotic region k* F r ^> 1 the first term xi( r ) 
dominates while X2 ( r ) is exponentially small. In particu- 
lar, it explains why the Xfi r ) term could not be obtained 
by the previous methodsa — the RKKY interaction as 
expressed in Ref.El was a series in powers of 1/r, and the 
exponential tail was evidently missing. 

On the contrary, the term X2(r) starts to play a de- 
cisive role at smaller distances k* F r < 1, where it pre- 



vails according to (27), (|29j). We stress that the con- 



dition k* F r > 1 was essential for the previous theories, 
while our expressions are applicable at weaker conditions 
(pp|) , (|2~l]) . Therefore the part X2( r ) represents the in- 
termediate asymptote of the RKKY interaction in the 
nearly-nested situation. 

Second we note that X2( r ) has the antiferromagnetic 
sign-reversal character. This feature of the RKKY inter- 
action at short distances appears to be-quite general. One 
can show for other types of dispersionllj that this kind of 
short range oscillations is always determined by the wave- 
vector connecting saddle points in the band structure. 

To clarify these two items further we consider the case 
of the perfect nesting of the electronic spectrum, /i = 
0, when the far asymptotic regime is not realized. We 
immediately see the disappearance of the first part of 
the RKKY interaction (|25|), while only this term could 
be obtained in the former methods of 1/r expansion. At 
the same time our second part of RKKY interaction ( l2§| ) 
survives. As a result we have 



X(r) 



1 



4ir 2 t \x 2 — y 2 



(30) 



This behavior indicates the log-squared singularity of the 
polarization operator on the antiferromagnetic wave vec- 
tor Q discussed, e.g., by Dzyaloshinskiilij. 

To validate our analytical findings (p5|), ( p§| ) we per- 
formed the direct numerical calculation of the RKKY 
interaction on the square \y\ < 10 with the spectrum 
given by (^). The results for /i = 0.1 and t — 0.5 are 
shown on the Figure We plotted on the Figure ||a 
the calculated value of the RKKY interaction versus the 
"distance" in the saddle point metrics, r* — y/\x 2 — y 2 \. 
According to this convention k* F r — {k* F ) max r* , so we 
can show the results for the whole plane in a simplest 
manner; for the chosen parameters {k* F ) max = \J /i/t ~ 
0.45. On the same plot we have drawn the curves 



xi(K 



X2{k* F r) and Xi{k* F r) + X2{k* F r), which are 



the predicted values of the interaction for the odd and 
even sites, respectively. No additional parameters were 
used. 



We see the remarkable agreement between the calcu- 
lated points and the theoretical formulas. As we expected 
at large distances the oscillations in the RKK Y are ob- 
served while at smaller distances r* < Jtj p, ~ 2.5 the 
situation is changed. The interaction for the "odd" sites 
x + y = 2n + 1 (i.e. for r = (0, 1), (1, 2), (0, 3). . . ) is 
of the antiferromagnetic (negative) sign. The interaction 
for the "even" sites (r = (0, 2), (1, 3), (0, 4) etc.) has a 
tendency to be ferromagnetic. In both cases the calcu- 
lated points closely follow our curves Xi( r ) ± X2( r ) up 
to the interatomic distances. The RKKY interaction ex- 
pected from the previously known expressions ( the term 
Xi(^F r )) is shown by a dashed line. 

To clearer represent the region of larger r* we mul- 
tiplied the calculated points x( r ) onto the appropriate 
values of (r*) 2 . The same was done for the theoretical 
curves, the results are shown on the Figure ||b. We see 
again that at large r* the interaction is characterized by 
the usually discussed asymptotic oscillations ( f27| ) . At the 
same time the difference between the "odd" and "even" 
sites is clear at the lower distances. 

Note that x( r ) f° r the diagonal x = ±y is not present 
on the Fig. @ and cannot be in principle compared to 
( p5| ) , ( |2"8| ) due to the criteria ( pp[ ) , ( |2l| ) ; we discuss it also 
in the next Section. 



B. noninteger values of R 

Let us now extend our analysis for the case of non- 
integer values of x, y. One can easily note that now 
the factors of the type exp[iir\x\sign(uj)} in (|22]), (|25|), 
( |A3| ) produce another source of discontinuity at ui = 0, 
in addition to the previously discussed one of the value 
ij2p(fx + iu>) . In particular both first and second terms 
in @ acquire the factors e ±2l7T ^ ^ 1 and e ±2l7T ^ ^ 1, 
respectively. Therefore both these terms now contribute 
although in a different manner. 

Consider first the case of \x\ > \y\ and \i < 0, which 
means p/i < and the clpsed character of the FS. A 
straightforward calculation^ shows then that the above 
expressions (p5|), (pgh are generalized as follows : 



Xi( r ) 



87rf 2 



[cos(27r|a;|)$i(fc^r) 



Xa(r) 



+ sin(27r|a;|)$^ 1) (fc>r) 
+ sm(27r\y\)^\k F r) 

ImI 



(31a) 



[cos(7r|a;| + n\y\)$ 2 (k* F r) 



AttH 2 

+ sin(7r|a:| + n\y\)^ A {k F r)} 



with the functions 



- \lY 2 (a) Jt(a)+Y?(a) J x 2 (a)] 
$i 2) (a) = ^[Kl{a)-Kl{a)] 



(31b) 

(32a) 
(32b) 
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$4 (a) = 



Yoia^ia) - Y 1 (a)K (a) 



(32c) G( n/2 ,n/2)(iu,r) 



The different terms appeared in (^T|) have different sig- 
nificance at large and small k* F r. 

In the far asymptotic regime k* F r > f the terms 



(2) 

$3 , $2, $4 are exponentially small and we find : 



X(r) cx 



sin(27r|a;| - 2k* F r) 



,r) 2 



> f 



(33) 



The period of oscillation in the above expression corr 
sponds to the notion of the calipering points on the FS 
We remind that these are the points where the direction 
of normal to the Fermi surface is (anti)parallel to the di- 
rection of r. In other words, the normal to the FS coin- 
cides with the direction of the Fermi velocity v = (v x , v y ) 
and it is parallel to r — r(cos tp, simp) provided 



u x sin (p — v y cos <p. 



(34) 



Near the saddle points (±7T, 0) one has v x /v y = 
—k x /ky and k 2 — k 2 — —fJ,/t. Therefore the caliper- 
ing points, satisfying the condition (|34|), are given by 
k c = ± (cos <p, — simp) [—(t/ n) cos2(^] -1 / 2 near the points 
(=F7r, 0), respectively. We measured the wave- vectors 
from the saddle points, therefore the true caliper of the 
fermi surface is given by the vector k c = (2tt, 0) — 2k c . 
The scalar product k c r is exactly what one finds in the 
Eq.(|3|) since k c r = ry/—ficos(2ip)/t = k* F r. 

At the smaller distances k* F r < 1 these are the terms 

<&z~\ < &3 2 ' ) , $2 7 which determine the main contribution to 
x(r). In this case one obtains : 



X(r) 



7rcos(7r|a;| + 7r|y|) + sin27r|a;| — sin 27r|y| 



47r 3 i|s 



(35) 



k* F r < I 



We see again that the interaction has a commensurate 
period of oscillations, although the oscillations for non- 
integer r are not described by a unique factor as it was 
in Eqs.@>), ©. 



V. THE FLAT PARTS OF THE SPECTRUM 

Let us discuss here the contribution to the RKKY 
interaction produced by the flat parts of the spectrum 
(|i"l|). First we observe that the formalism developed in 
the main part of the paper cannot be applied to the 
vicinities of the points (±7r/2, ±7t/2) since all the com- 
ponents of the mass tensor are infinite at these points. 
This very special case should be treated separately; one 
can also distinguish here the regions of intermediate and 
far asymptotes. At the intermediate distances r < \t/ 
we obtain the Green's function from the vicinities of 
(7r/2,7r/2) and (-tt/2, -tt/2) in the form 



2t 



-8 K (x — y) exp [i|x|(7r + z /2t)sign(uj)] 

(36) 



and the corresponding Green's function from the vicini- 
ties of (7r/2, —tt/2) and (—tt/2, tt/2) is obtained from this 
expression by the replacement y — > —y. The function 
5 K (x) in ( |36| ) has the 5— function-like properties and is 
defined by 



siiiki 1 

K (X) = , K ~ - 

ttx a 



(37) 



We wish to point out that the power-law decrease of S K (x) 
at large x stems from our assumption that the absence 
of dispersion along x is lost abruptly at \k x \ > k. In 
fact the expression ( |37| ) is the Fourier transform of 6(n — 
\k x \). In general the dependence of dispersion on k x is 
much smoother ; as a result, the decay of 5 K (x) at large x 
should be much faster, while the 5— function-like property 
preserves. 

We see that the above Green's function has a sizeable 



< f. Outside this domain 
(0) is 



values only in a band \x — y\ 
the principal contribution to the total G(iu,r) 
delivered by G( i7r )(iw, r) and G^^iiu), r), Eqs. < \22\ ) and 
|). Particularly it means (see Fig.|l|) that the terms of 



(I 

the type G^/2^/2) G(o,tt) m the expression ( ^4[ ) should 
not be considered. 

As a result the contributions to the RKKY interaction 
from the flat parts of FS acquire the following form : 

XfUr) = ^"tll^ W(* -v) + + y)} - (38) 



Airt\x\ 

< I*/mI 



here two terms in the square brackets correspond to the 
different regions in the r-space. We see that this part of 
interaction which is present along the diagonals is slo wly 
decaying as I/r. The amplitude of it, according to ( |37| ) 
has the model cutoff parameter k 2 . Hence we cannot di- 
rectly compare this part of RKKY with the results of our 
numerical calculations, although the overall 1/r depen- 
dence of RKKY interaction along the diagonals and slow 
oscillations are verified by the numerical data as well. 
At small integer values of x — y the RKKY term fl3S| ) 
corresponds to the ferromagnetic sign of the interaction 
between the localized moments. This behavior however 
does not define the particular type of magnetic ordering 
and it is the term ( |30| ) which determines it. 

At extremely large distances close to diagonal r ^ 
\t/fi\ the dropped /c 3 — terms in the expansions (|l0|), ( fTl] ) 
become important. Hence the spectrum becomes essen- 
tially two-dimensional, with the corresponding change in 
the character of RKKY. Near the diagonals ip = ±7r/4+<j!> 
one has : 



Xfiat(r) = - 



sin|a;|(27r + ^A(l - (j) 2 /^)) 
A7r 2 \^\x 2 [l + ^ 2 / ( j ) 2 ] 



t 

Ml 
(39) 



G 



This far asymptote of RKKY interaction from the flat 
parts of dispersion holds in the narrow sectors near the 
diagonals \(j>\ < (j> = V2\n/U\, ( cf. ©, © and Fig. 
[j] ). It has the 1/r 2 dependence while its period of os- 
cillations correponds to the notion of calipering points 
discussed above. 
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VI. CONCLUDING REMARKS 



APPENDIX A: THE STATIONARY POINT AT 
THE ZONE BOUNDARY 



It is worthwhile to compare our expressions for the 
Green's .function with previous results. It was observed 
(see e.g.EJ) that for the tight-binding spectrum (0) one 
can use the recurrence relations to express the value of 
G(v,r) in terms of G(w,0). It was noted also however 
that using these relations one meets the numerical in- 
stabilities at large r. Alternatively G(u>, r) can be esti- 
mated at large r by the steepest descent method. The 
solution obtained by this latter method corresponds to 
the asymptotes of Eqs. fll7j), (|l8|) for the case of large 



In this sense our expressions extend the 
previous findings for the Green's function and provide 
the analytical formulas in the region of the intermediate 
distances 1 < r < l/k* F . 

Let us briefly discuss here the role of finite temper- 
atures for our treatment. In this case instead of the 
integral ( pi] ) one considers the sum over the Matsub- 
ara frequencies with the Green's functions given by 
(p2|), (p3|). With the use of analytical continuation, 
this sum can be represented as the integral over the 
real axis of lu. One can note however directly from 
the form of the Green's functions that the effect of fi- 
nite temperatures is important when T exceeds the ef- 
fective chemical potential p. At large distances r > £ = 

( k *F )~ 1 Vv/ T K V / V^ W&ihave the RKKY interaction 
exponentially suppressed .L3 The opposite case r % £ cor- 
responds essentially to the case p = described by the 
Eq. (|30|) . Therefore the far asymptote of RKKY leading 
to possible incommensurate magnetic ordering is absent 
in this case and we remain with the only tendency to 
commensurate AF order. 

In conclusion we found the closed analytic expres- 
sions for the RKKY interaction in a layered metal with 
nearly nested Fermi surface. Along with the usual 2k F - 
oscillations realized at far distances we demonstrate the 
existence of the intermediate asymptote of the interac- 
tion. This latter asymptote has the commensurate AF 
type of oscillations and is the only term surviving at the 
exact nesting. We show that our analytical formulas are 
in the good accordance with the numerically found values 
of interaction in a range up to near interatomic distances. 
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Let us consider the case of the stationary point lying on 
the boundary of the Brillouin zone. We discuss it on the 
example of the sum of the points (tt, 0) and (— tt, 0) for the 
tight-binding form of dispersion (j?]). The integration over 
k y meets no difficulties, while the domains of integration 
over k x , namely (— tt, — 7r + k) and (tt — k,tt), could be 
combined as follows : 



r K ^ 

G Z B(iu,x) = e~ mx \ — -exp 
Jo 27T 

+ (x — > —x) 
= A 



k 2 



2to 7 



. k*y x 

cos(7rx)Eri i— 

2 7 



+2i sin(7ra) Erf ( —i 



with 



A = — exp 

7T7 



2t 



(Al) 



(A2) 



2t 



e*" — , ReOy) > 0. 



The demand for the argument ±^7/2 — ix/"f in (Al) 
to be large corresponds to the criteria of applicability 
of our expressions (^fj), (|2l]) and we do not discuss it 
here. Similar to the above treatment at k* F r 3> 1 one has 
t ~ k* F r/(2n) and \x/j\ > k* F r > 1. At k* F r < 1 we 
saw that \p\ < r < \z\ ; on the other hand the values of 
\z\ at which the integral (|I|) saturated were of order of 
l/\p\. Therefore we can let r ~ \p\ and obtain \x/j\ ~ 
yj\x 2 m x / p\ > 1. Hence in b oth cases the argument of 
the second error function in (Al) is large enough, while 
its sign coincides with the sign of a product —xojm x . 
Therefore the Green's function has a following form : 



G z b(iuj,x) — v / 7r^4exp[— iir\x\sign(umi x )] 



(A3) 



This result means that the Gaussian value of integral 
is attained in one of the above domains of integration. 
Note that G z b(iuJ,x) is an even function of x ; it is a 
manifestation of the reflecting property of the Brillouin 
zone boundaryO. For a particular choice of our spectrum 
(||) we have m x = —l/2t < and thus the formula (23). 
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FIG. 1. The regions of applicability of the different 
Green's functions in the sector x > 0, y > of the r-space. 
The contributions from the saddle points of the spectrum are 
given by (|22|, (0) and are defined away from the diagonal. 
The expression (B6|) provides the contribution from the flat 
parts of dispersion. The solid line is the border between the 
applicability regions of corresponding equations. At the same 
plot the curve k* F r — 1 is shown by the dashed line. 



* 



0.00 



u 

K -0.05 





(b) 


\ 4 \ 

\ tr V\ 




\ o f v\ 
\ ,/ X\ 

\ M v\ 




,1 w 

97 \ 
x t / 

v -- y 1 




J odd sites 




o o / • even sites 





1 2 3 4 5 6 7 

1/2 



r =|x2-y2 |1/ 



8 



FIG. 2. a) The calculated RKKY interaction for the 
tight-binding spectrum ((?]) with fi — 0.1, t = 0.5 is shown 
along with the theoretically obtained curves, b) The same 
quantity plotted in a modified way, for clearer representa- 
tion of the region of larger r*. The actual coordinates of 
some sites are indicated. One can note the remarkable ac- 
cordance between the numerical findings and the analytical 
results in the intermediate region, where the curves for "odd" 
and "even" sites differ visibly. The asymptotic solution for 
RKKY extended to shorter distances is shown by a dashed 
line for comparison. 
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